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ABSTRACT 

For early-type galaxies, the correlations between stellar mass and size, velocity dis- 
persion, surface brightness, color, axis ratio and color-gradient all indicate that two 
mass scales, M* — 3 x 10 10 M Q and M* = 2 x 10 11 M Q , are special. The smaller scale 
could mark the transition between wet and dry mergers, or it could be related to the 
interplay between SN and AGN feedback, although quantitative measures of this tran- 
sition may be affected by morphological contamination. At the more massive scale, 
mean axis ratios and color gradients are maximal, and above it, the colors are redder, 
the sizes larger and the velocity dispersions smaller than expected based on the scaling 
at lower M*. In contrast, the color-a relation, and indeed, most scaling relations with 
a, are not curved: they are well-described by a single power law, or in some cases, are 
almost completely flat. When major dry mergers change masses, sizes, axis ratios and 
color gradients, they are expected to change the colors or velocity dispersions much 
less. Therefore, the fact that scaling relations at a > 150 km s~ show no features, 
whereas the size-M*, 6/a-M*, color-M* and color gradient-M* relations do, suggests 
that M* = 2 x lO n M0 is the scale above which major dry mergers dominate the 
assembly histories of early-type galaxies. 
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1 INTRODUCTION 

Recent work (Bernardi et al. 2010b) has shown that the 
color-magnitude relation of early-type galaxies in the SDSS 
differs significantly from a pure power-law, curving down- 
wards at low and upwards at large luminosities (A/ r > —20.5 
and M r < —22.5, respectively). This is also true of the color- 
size relation, and is even more apparent with stellar mass, 
where the corresponding mass scales are M, < 3 x 10 10 M o 
and M* > 2 x 10 Mq, respectively. The upwards curva- 
ture at the massive end does not appear to be due to stel- 
lar population effects. Curvature in the color-luminosity (or 
color- M*) relation at the faint end was noticed before (e.g. 
Graham 2008; Skelton et al. 2009); the curvature at bright 
end is the new finding of Bernardi et al. (2010b). 

Curvature at the low mass end using other parameters 
(e.g. surface brightness) has been noticed before as well (e.g. 
Kauffmann et al. 2003; Shankar et al. 2006); the subject of 
this Letter is to analyse non-linear scaling relations at the 
high mass end, extending the analysis of Hyde & Bernardi 
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(2009). We show that the curvature in the color— M, re- 
lation coincides with curvature in other relations with M*. 
However, when M* is replaced with velocity dispersion, then 
there is little curvature. Although curvature in scaling re- 
lations does not imply a change in the physics which sets 
the relations (e.g. Graham & Guzman 2003; Graham 2010), 
we argue that our findings suggest that major dry mergers 
dominate the mass growth at M* > 2 x 10 n MQ. 

Our results are based on two different ways of select- 
ing early-type samples from the SDSS database. One fol- 
lows Hyde & Bernardi (2009): the image must be round 
(b/a > 0.6) and the light profile shape must be well-fit by 
a deVaucoleurs profile (fracDev = 1). The other is a sim- 
ple cut on how centrally concentrated the surface bright- 
ness is (C r > 2.86). The former method produces a sam- 
ple that is more purely elliptical; the latter contains many 
edge-on disks. See Bernardi et al. (2010a) for a more de- 
tailed discussion of these selection criteria, and of the SDSS 
photometric and spectroscopic parameters which we use 
below. Where necessary, we have assumed a spatially flat 
background cosmology with energy density dominated by 
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Figure 1. Curvature in the correlations between stellar mass and (from top to bottom, left) size, velocity dispersion and surface 
brightness, and color, axis ratio, and color gradient (top to bottom, right), in the Hyde-Bernardi sample. The vertical dashed lines mark 
the scales where some of the relations change slope: M* = 3 X 10 10 Mq and 2 X 10 11 Mq, which correspond approximately to M r = —20.5 
and —22.5. Scalings in a sample selected to have C r > 2.86 are shown only where they differ from the scalings in the Hyde-Bernardi 
sample. 



a cosmological constant A = 0.7, with a Hubble constant 
Hq = 70 km s _1 Mpc -1 at the present time. 



2 CURVATURE IN RELATIONS WITH M», 
BUT NOT WITH a 

Figure [T] shows correlations between stellar mass and (from 
top to bottom, left) size, velocity dispersion and surface 
brightness, and color, axis ratio, and color gradient (top 
to bottom, right), in the Hyde-Bernardi sample. (We dis- 
cuss how we define the gradient in Section 13.21 ) None of 
these correlations are pure power laws. Although the curva- 
ture at log 10 M*/Mq < 10.5 is interesting - galaxies at the 
faint, low mass end (M r > -20.5, log 10 (M*/M Q ) < 10.5) 
tend to curve towards bluer colors, larger sizes, fainter sur- 
face brightnesses, smaller axis ratios and color gradients - 
in what follows, we will focus on what appears to be a tran- 
sition mass scale at higher masses. At log 10 M*/Mq > 11.3, 



the relations curve towards larger sizes, smaller than ex- 
pected velocity dispersions, fainter surface brightnesses (left 
panels), and smaller axis ratios and smaller color-gradients 
(right panels). This is precisely the mass scale on which the 
color-M* relation curves towards redder colors (top right 
panel). 

Figure [2] shows that when Af» is replaced with velocity 
dispersion, then there is little curvature at logcr/kms -1 > 
2.2. In fact, the correlations with surface brightness, color 
gradient and axis ratio are almost completely flat. (That 
surface brightness and a are uncorrelated was noted by 
Bernardi et al. 2003.) The fact that there is no feature at 
the largest a in any of these relations, despite clear features 
in the scalings with M», is what has motivated this Letter. 

In this context, it is important to note that the relation 
between Md yn oc Ra 2 and luminosity or M* is very well 
described by a single power-law over the entire range: the 
curvature in the sizes and velocity dispersions cancel (Fig- 
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Figure 2. Curvature in the correlations between velocity dispersion and (from top to bottom, left) size, stellar mass and surface 
brightness, and color, axis ratio, and color gradient (top to bottom, right), in the Hyde-Bernardi sample. The dashed lines mark the 
scales where one would expect to see a change in the slope of the relations based on Figure [T] Scalings in a sample selected to have 
C r > 2.86 are shown only where they differ from the scalings in the Hyde-Bernardi sample. 



ure[3j)- Presumably, this is because the objects we observe 
are virialized, whatever their merger histories. 



for axis-ratios and color-gradients, and, significantly, that 
the curvature is absent when M* is replaced by a. 



3 DISCUSSION 

Major dissipationless mergers are expected to change the 
sizes in proportion to the masses, but to leave the velocity 
dispersions and colors unchanged. In contrast, minor dissi- 
pationless mergers produce larger fractional changes in size 
than in mass, and decrease the velocity dispersions and col- 
ors (see Appendix C in Bernardi et al. 2010b for details). 
Therefore, the curvature in the correlations between M* and 
other parameters as size, a and color, which have been no- 
ticed before, have all been discussed in this context (e.g. 
Davies et al. 1983; Matkovic & Guzman 2005; Bernardi et 
al. 2007; Hyde & Bernardi 2009; Bernardi et al. 2010a,b). 
What is new here is the recognition that these all occur at 
the same mass scale, that this mass scale is also important 



3.1 Axis-ratios 

The b/a — M* relation (right center panel of Figure [1} 
deserves further comment. Van der Wei et al. (2009) re- 
port that the width of the b/a distribution changes at 
log(M*/MQ) ~ 10.5. They interpret this as evidence that, 
above this mass, assembly histories are dominated by major 
mergers. Our results suggest this is not the full story. 

In Figure [1] we have shown two versions of this relation, 
because the Hyde-Bernardi selection requires b/a > 0.6. In 
this sample, b/a decreases at log(M*/M0) > 11.3. However, 
notice that this decrease is even more marked in the sample 
selected to have CY > 2.86, where no cut on b/a is applied. 
Compared to the Hyde-Bernardi sample, this sample has 
considerably smaller b/a at small Af, . Bernardi et al. (2010a) 
show that this is primarily due to an increased incidence 
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of disks and contamination by Sas, because the CV > 2.86 
sample is not as purely elliptical/early- type as the Hyde- 
Bernardi sample. We believe this change in morphological 
mix is the primary reason why Van der Wei et al. saw what 
they did. 

We believe the real feature of interest is the drop in b/a 
at log(M»/M Q ) > 11.3 where (Bernardi et al. 2010a show 
that) morphological mix is no longer an issue. Van der Wei et 
al. also see this drop, but they dismiss it. Instead, we believe 
the narrowing of the distribution at log(M*/Mo) ~ 10.5 
marks the transition from dissipational to dissipationless 
histories, or a change in relative importance of SN and AGN 
feedback (e.g. Kauffmann et al. 2003; Shankar et al. 2006), 
while the decrease in b/a at log(Af*/M0) > 11.3 marks the 
transition to major dry mergers. This decrease has been ex- 
pected for some time (see Gonzalez-Garci'a & van Albada 
2005; Boylan-Kolchin et al. 2006; Ragone-Figueroa & Plio- 
nis 2007; Ragone-Figueroa et al. 2010) - it was first found 
by Bernardi et al. (2008). This is thought to indicate an in- 
creasing incidence of major radial mergers, since these would 
tend to result in more prolate objects. 

3.2 Color gradients 

The right bottom panel of Figure \T\ shows that color- 
gradients — here defined to be the difference between 
the model and Petrosian colors - - are maximal at 
log(M»/M0) ~ 11.3. (The former are approximately the 
color within the half-light radius, whereas the latter are more 
closely related to the ratio of the total luminosities in g and 
r, so correspond to larger scales.) This is consistent with 
Figure Bl of Bernardi et al. (2010b), who used the same 
definition of color gradient. It is also consistent with Roche 
et al. (2010), who used a different estimator of the gradient: 
the ratio of the half-light sizes in the g and r bands. As we 
discuss below, we believe that the appearance of this same 
mass scale is again signaling the onset of major dry mergers. 

Whereas major mergers are expected to decrease color 
gradients (e.g. Di Matteo et al. 2009), minor mergers should 
not change the gradients significantly (Kobayashi 2004) or 
they may enhance them slightly. This is because the smaller 
bluer object involved in the minor merger is expected to 
deposit most of its stars at larger distances from the center of 
the object onto which it merged. (If it had its own gradient, 
then the bluest of its stars would have been deposited at the 
largest radii.) 

As a simple check of this argument, note that major 
mergers, which double M, , do not change a (Appendix C 
of Bernardi et al. 2010b). Therefore, a plot of color gradient 
versus a should show less of a feature than when gradients 
are plotted versus M„ . This is indeed what we see in Figuref2] 
(right bottom panel) - the correlation is almost flat at a > 
150 km s . Such a plot should also show greater scatter, 
since a range of merger histories, hence gradients, can all 
have the same a. While we do see this increase in scatter, 
note that the scatter in the gradient- M» relation grows even 
more dramatically - something that is not easily explained. 
On the other hand, the major merger picture also provides 
a natural explanation for why none of the scaling relations 
in Figure [2] show any feature at logcr/kms -1 > 2.2, and 
those that are most clearly sensitive to merger histories are 
almost completely flat. 
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Figure 3. No curvature in the correlations between dynamical 
mass and luminosity (top) and stellar mass (bottom). 

Mergers are not the only way to produce or alter color 
gradients. In some models, gradients are related to feed- 
back and winds (Pipino et al. 2010). Our demonstration 
that gradients scale differently with M* than with a may 
have interesting implications for such models. In addition, 
producing the downturn we see at log(M*/M0) > 11.3 is an 
interesting challenge for such models, as is relating this to 
the changes in size, a and b/a we have found. Because the 
major dry merger model provides a simple framework for 
understanding all these relations, our results suggest that 
M t > 2 x 1O 11 M0 is the scale above which major dry merg- 
ers dominate the assembly history. 



4 IMPLICATIONS 

We have found that a variety of early-type galaxy scaling re- 
lations - the size-M*, b/a-M t , color- M t and color gradient- 
M* relations - all show departures from a pure power-law at 
M t — 2 x 1O 11 M0 , whereas there is no such feature when M» 
is replaced with a. Since major dry mergers are expected to 
change the sizes, axis ratios and color gradients of galaxies 
while leaving the velocity dispersion and color unchanged, 
our findings suggest that the total stellar mass in early-types 
with M* > 2 x lO 11 Af0 today must have grown primarily 
by relatively recent major dry mergers. In Bernardi et al. 
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(2010b), we argue that such mergers may be required to 
reconcile the z ~ 1 counts of objects with M 4 > 2 x 10 11 Mq 
with those at z ~ 0. 

This particular mass scale also appears in analyses of a 
local sample (higher quality data, but significantly smaller 
sample), where it is identified with the transition to dry 
mergers (see page 270 and related discussion in Kormendy 
et al. 2009). It is special in hierarchical models also. Figure 
3 of Guo & White (2008) shows that below 1.6 x lO n M 
star-formation has been a significant part of the mean stellar 
mass growth rate (much of it through wet mergers), whereas 
the stellar mass growth at masses above this occurs only 
through dry mergers. See Hopkins et al. (2008) and Eliche- 
Moral et al. (2010a, b) for other arguments suggesting dry 
mergers since z ~ 1 are a natural and necessary part of the 
assembly history at M* > 2 x 10 11 Mq. While it is reassuring 
that many lines of study all identify this same mass scale, 
we feel it worth emphasizing that our analysis suggests that 
above this mass scale, the mergers were not just dry - they 
were major. In this respect, there is some tension between 
our conclusions, and recent work which argues that although 
the mass in the central kpc or so of early-type galaxies has 
not grown since z ~ 2, the half-light radii have increased by 
more than a factor of two. This suggests that, since z < 2, 
mass has been added to the outer regions only: this sort of 
inside-out scenario for the growth is most easily understood 
if the mergers were minor (e.g. Lapi & Cavaliere 2009; Cook 
et al. 2009; Bezanson et al. 2009). However, as Tiret et al. 
(2010) note, the observation of constant mass in the central 
regions does not, by itself, exclude major mergers. In the 
simulations of Gao et al. (2004), as an object assembles its 
mass through major dry mergers, the mix of particles in 
the central regions can change dramatically, even though 
the total mass in the central regions remains constant. Our 
finding that color-gradients are erased at large masses may 
be indicating that this is indeed what happens at M* > 
2 x 10 n M Q . 

Finally, it is interesting to ask how BCGs, which are 
amongst the most massive objects in the local universe, fit 
into this picture? Compared to non-BCGs of similar mass 
or luminosity, their colors are slightly redder (Roche et al. 
2010; Figure 10 in Bernardi et al. 2010b), they have smaller 
color gradients (Roche et al. 2010), and slightly larger sizes 
(Bernardi 2009). Whereas the first two are in agreement with 
our major merger picture, the last suggests more size growth 
than is usually associated with major mergers. Hence, it may 
be better to think of BCG formation as a two step process. 
In the first, the major mergers which result in the object 
becoming a BCG erase its color gradient (and decrease b/a 
- Bernardi et al. 2008); thereafter, minor mergers puffed up 
its size. Tidal stripping during the minor merger may also 
have contributed to the formation of intracluster light in 
its host halo (Bernardi 2009; Bernardi et al. 2010b). This 
two step picture is in striking agreement with a detailed 
analysis of the age, metallicity and abundance gradients of 
BCG NGC 4889 (Coccato et al. 2010). 
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